Differential cross sections and spin density matrix elements for the reaction *yp — ► puo 
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High-statistics differential cross sections and spin density matrix elements for the reaction 7p — > pui 
have been measured using the CEBAF large acceptance spectrometer (CLAS) at Jefferson Lab for 
center-of-mass (cm.) energies from threshold up to 2.84 GeV. Results are reported in 112 10- 
MeV wide cm. energy bins, each subdivided into cos#" m . bins of width 0.1. These are the most 
precise and extensive u photoproduction measurements to date. A number of prominent structures 
are clearly present in the data. Many of these have not previously been observed due to limited 
statistics in earlier measurements. 

PACS numbers: 11.80.Cr,11.80.Et,13.30.Eg,14.20.Gk,25.20.Lj,23.75.Dw 



I. INTRODUCTION 

Studying low-energy lo photoproduction presents an 
interesting opportunity to search for new baryon res- 
onances. Previous experiments have produced cross- 
section measurements with relatively high precision at 
most production angles; however, precise spin density 
matrix elements have only been measured at very forward 
angles [TJ |3J HI [5] . Theoretical interpretation of these 
data indicate strong i-channel contributions from both 7r° 
and Pomeron exchange, while the backward peak in the 
cross section has been interpreted as evidence of nucleon 
u-channel contributions [5J [51 [!]• Several attempts 
have been made to extract resonant contributions that 
have obtained conflicting results [8| ITO l ITT ] [T2 | [T3 | H4], 
Precise polarization information is needed in order to 
place stringent constraints on the physics interpretation 
of lo photoproduction data. 

The impact of polarization information can be seen 
by comparing the partial wave analysis results ob- 
tained using only cross-section data [13] to those that 
also included the low-precision polarization results from 
SAPHIR [H]. The former found that at threshold the 
dominant s-channel contributions are from the Pi3(1720) 
and Fl 5 (1680), while the latter found that the £>i 5 (1675) 
and i 7 'i5(1680) are dominant in this region. Including po- 
larization information, even with very limited precision, 
provided strong additional constraints on the interpreta- 
tion of the data. Thus, obtaining high-precision polariza- 
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tion results is a vital step towards understanding baryon 
resonance contributions to lo photoproduction. 

Beyond this, quark model calculations of baryon de- 
cays |15j predict that a number of the so-called missing 
baryons should couple to coN final states. In particu- 
lar, in the above model, nearly all of the missing posi- 
tive parity N* states are expected to have non-negligible 
couplings to loN. Thus, good data on lo photoproduction 
coupled with a partial wave analysis (PWA) could pro- 
vide important new information on light-quark baryons. 

The data presented here are part of a larger program 
to simultaneously measure photoproduction of mesons 
and then carry out partial wave analyses on the resulting 
data. This article presents differential cross section and 
spin density matrix element measurements for lo photo- 
production. In a companion article published concurrent 
to this one [16j, we present a detailed partial wave analy- 
sis of these data where clear s-channcl resonance contri- 
butions are identified. A forthcoming article will discuss 
the impact on current theoretical models and coupled- 
channel analyses of these new precise measurements |17j . 



II. EXPERIMENTAL SETUP 

The data were obtained using the CEBAF large ac- 
ceptance spectrometer (CLAS) housed in Hall B at the 
Thomas Jefferson National Accelerator Facility. Real 
photons were produced via bremsstrahlung from a 4 GeV 
electron beam hitting a 10~ 4 radiation length gold foil. 
The recoiling electrons were then analyzed using a dipole 
magnet and scintillator hodoscopes in order to obtain, or 
"tag" , the energy of the photons [IB] (the so-called pho- 
ton tagger). The tagging range and energy resolution 
were 20% — 95% and 0.1% of the electron beam energy, 
respectively. The useful center-of-mass (cm.) energy 
(W) range for this analysis was from w-photoproduction 
threshold at W = 1.72 GeV up to 2.84 GeV. In this 
range, the data were analyzed in 10-MeV wide W bins. 

The physics target, which was filled with liquid hydro- 
gen, was a 40-cm long cylinder with a radius of 2 cm. 
Continuous monitoring of the temperature and pressure 
permitted determination of the density with uncertainty 
of 0.2%. The target cell was surrounded by 24 "start 
counter" scintillators that were used in the event trigger. 
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FIG. 1: (a) The confidence levels resulting from four- 
constraint kinematic fits performed on a sample of events to 
the hypothesis ~/p — > p7r + 7r~ integrated over all kinematics. 
The "peak" near zero consists of events that do not match 
the hypothesis, along with poorly measured (due to multi- 
ple scattering, etc.) signal events. Agreement with the ideal 
(flat) distribution for signal events is very good, (b) Example 
pull-distribution for the momentum of the 7r~ from the same 
kinematic fits as in (a). Only events with a confidence level 
larger than 1% are shown. The line represents a Gaussian 
fit to this distribution. For this event sample, the parame- 
ters obtained are fi = -0.029 ± 0.001, a = 1.086 ± 0.001 (the 
uncertainties are purely statistical), which are in very good 
agreement with the ideal values /j, = 0, a = 1. Both (a) and 
(b) are good indicators that the CLAS error matrix is well 
understood. 



The CLAS detector utilized a non-uniform toroidal 
magnetic field of peak strength near 1.8 T in conjunc- 
tion with drift chamber tracking to determine particle 
momenta. The detector was divided into 6 sectors, such 
that when viewed along the beam line it was six-fold sym- 
metric. Charged particles with laboratory polar angles in 
the range 8° — 140° could be tracked over approximately 
83% of the azimuthal angle. A set of 288 scintillators 
placed outside of the magnetic field region were used in 
the event trigger and during offline analysis in order to 
determine time of flight (TOF). The momentum resolu- 
tion of the detector was, on average, about 0.5%. Other 



components of the CLAS, such as the Cerenkov counters 
and the electromagnetic calorimeters, were not used in 
this analysis. A more detailed description of the CLAS 
can be found in Ref. [19] . 

The event trigger required a coincidence between sig- 
nals from the photon tagger and the CLAS. The signal 
from the tagger consisted of an OR of the first 40 of 
the 61 total timing scintillators, corresponding to photon 
energies above 1.5 GeV. Recording of events associated 
with photons hitting counters 41-61 required a random 
tagger hit in counters 1-40 during the trigger timing win- 
dow. This allowed for the acquisition of greater statistics 
at the higher photon-energy range of this experiment. 
The signal from the CLAS required at least two sector- 
based signals. These signals consisted of an OR of any 
of the 4 start counter scintillators in coincidence with an 
OR of any of the 48 time-of-flight scintillators in the sec- 
tor. The rate at which hadronic events were accumulated 
was about 5 kHz; however, only a small fraction of these 
events contained the reaction of interest to the analysis 
presented here. 



III. DATA AND EVENT SELECTION 

The data used were obtained in the summer of 2004 
during the CLAS "glla" data taking period, in which ap- 
proximately 20 billion triggers were recorded. The rela- 
tively loose electronic trigger led to accumulation of data 
for a number of photoproduction reactions. The relative 
timing of the photon tagger, the start counter and the 
time-of-flight elements were aligned during offline calibra- 
tion. Calibrations were also made for the drift times of 
each of the drift chamber packages and the pulse heights 
of each of the time-of-flight counters. Finally, processing 
of the raw data was performed in order to reconstruct 
tracks in the drift chambers and match them with time- 
of-flight counter hits. 

The reconstructed tracks were corrected for small im- 
perfections in the magnetic field map and drift cham- 
ber alignment, along with their mean energy losses as 
they passed through the target, the beam pipe, the start 
counter and air. In addition, small corrections were made 
to the incident photon energies to account for mechanical 
sag in the tagger hodoscope. 

The CLAS was optimized for detection of charged par- 
ticles; thus, the 7T + 7r~7r° decay of the u> was used to select 
the reaction of interest in this analysis. Detection of two 
positively charged particles and one negatively charged 
particle was required. A cut was placed on the confi- 
dence levels obtained from one-constraint kinematic fits 
performed to the hypothesis jp — > pir + ir~ (ir°) in order to 
select events consistent with a missing ir° . All negatively 
charged tracks were assigned a 7r~ identification. Kine- 
matic fits were run for each of the possible p, tt + particle 
assignments for the positively charged tracks using each 
of the recorded photons in the event. The combinations 
of photons and charged particles with confidence levels 
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greater than 10% were retained for further analysis. 

The covariance matrix was studied using four- 
constraint kinematic fits (energy and momentum con- 
servation imposed) performed on the exclusive reaction 
7P — * pTT + 7r~ in both real and Monte Carlo data samples. 
The confidence levels in all kinematic regions were found 
to be sufficiently flat and the pull-distributions (stretch 
functions) were Gaussians centered at zero with a = 1 
(see Fig. [TJ. The uncertainty in the extracted yields due 
to differences in signal lost because of this confidence- 
level cut in real as compared to Monte Carlo data is es- 
timated to be 3% -4%. 

The tagger signal time, which was synchronized with 
the accelerator radio-frequency (RF) timing, was propa- 
gated to the reaction vertex in order to obtain the start 
time for the event. The stop time for each track was ob- 
tained from the TOF scintillator element hit by the track. 
The difference between these two times was the measured 
time of flight, t meas . Track reconstruction through the 
CLAS magnetic field yielded both the momentum, p, of 
each track, along with the path length, L, from the reac- 
tion vertex to the time-of-flight counter hit by the track. 
The expected time of flight for a mass hypothesis, m, is 
then given by 



passing the fit under the proton hypothesis. The re- 
gion near (0, 0) contains events where both tracks are 
good matches to their respective particle identification 
hypotheses. The 2-ns radio-frequency time structure of 
the accelerator is evident in the out-of-time event clus- 
ters. Events outside of the black lines, where neither hy- 
pothesis was met, were cut from our analysis. This cut 
was designed to remove a minimal amount of good events. 
The Feldman-Cousins method EQ] was used to place an 
upper limit on the signal lost at 1.3%. Any remaining 
accidental events fell into the broad background under 
the uj, and were rejected during the signal-background 



separation stage of the analysis discussed in Sec. IV 



u exp 




(1) 



Fiducial cuts were applied on the momenta and an- 
gles of the tracks in order to select events from the well- 
understood regions of the detector. Included in these 
cuts was the removal of 13 of the 288 time-of-flight el- 
ements due to poor performance. In addition, events 
where the missing ir° was moving along the beam line, 
cos > 0-99; were cut in order to remove leakage from 
the jp — > ptt + tt~ reaction. A more detailed description 
of the analysis procedures presented in this article can be 
found in Ref. El]. 



IV. SIGNAL-BACKGROUND SEPARATION 



In addition to uj events, the resulting sample consisted 



The difference in these two time-of-flight calculations 

Atof = t meas - t exp , was used in order to separate pro- of events from thc reaction 7 p 
tons from pions and to remove events associated with 
out-of-time photons. 




A tof proton (ns) 



FIG. 2: (Color Online) Aio/„.+ (ns) versus Ato/ P (ns): Parti- 
cle identification cut for a sample of events that pass a 10% 
confidence level cut when kinematically fit to the hypothesis 
IP pir + ir~ (ir°). The black lines indicate the timing cuts. 
Note the logarithmic scale on the intensity axis. 

Fig. [2] shows Atof for the track passing the kinematic 
fit under the 7r + hypothesis versus Atof for the track 



where the 

three-pion invariant mass was consistent with that of the 
uj. These background events could arise from a variety 
of physics reactions, but they all share the characteristic 
that they cannot reproduce the narrow mass structure 
associated with the uj. Typical background levels were 5 
to 10% of the uj peak height, while in a small number of 
bins near threshold and where the cross section is very 
low (e.g. near cos6^ m of —0.2 at the highest photon 
energies), the background was as large as 25% of the u> 
peak. Thus, the three-pion invariant mass could be used 
as a tool to separate the signal from the background. 

In order to accurately extract physical observables for 
uj photoproduction, background events (all non-w events) 
must be separated from the signal in a way that pre- 
serves all kinematic correlations. The method we have 
developed, described in detail elsewhere ED E2], was 
used to assign each event a signal weight factor, Q, or 
equivalcntly, a background weight factor, 1 — Q. These 
Q- factors were then used to weight each event's contri- 
bution to the "log likelihood" during the event-based fits 
discussed in Sec. [V] These fits were used to obtain the 
detector acceptance and to extract the spin density ma- 
trix elements. The Q- factors were also used to weight 
each event's contribution to the differential cross section. 

The key feature of our procedure involves selecting 
each event's N c "nearest neighbor" events (we chose 
N c = 500). This is done by defining a metric in the 
space of all relevant kinematic variables, excluding the 
37r invariant mass, M^. Each subset of events occupies 
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FIG. 3: (Color Online) (a) The 3-7T invariant mass distri- 
bution in the W — 2.205 GeV bin, integrated over all kine- 
matics, for all events (unshaded) and for events weighted by 
the background factors, 1 — Q (shaded), (b) The A distribu- 
tion of events in the same W bin that satisfy {M^^ — M u \ < 
25 MeV/c 2 (unshaded), the same events weighted by signal 
factors Q (dashed-red) , and by background factors 1 — Q 
(shaded) . The line represents a fit of the signal to the function 
a\. 



a very small region of phase space; thus, the M 3v dis- 
tribution can safely be used to determine each event's 
Q-f actor, while preserving the correlations present in the 
remaining kinematic variables. This method facilitates 
separation of the signal and background without having 
to resort to dividing the data up into bins. Binning the 
data is undesirable due to the high dimensionality of the 
reaction being studied in this analysis. 

To this end, unbinned maximum likelihood fits were 
carried out for each event, using its nearest neighbors, 



to determine the parameters a = (bo, b±, &2, b 3l b^, s, a) in 
the probability density function 



F(M 37ri a) 



B(M 3 *,a) + S(M 3v ,a) 
J (B(M 3 „,a) + S(M 3 „,a)) dM 3n 



where 



S(M 3 „ , 5) = s ■ V (M 3n , M u , T w , a) 



, (2) 



(3) 



parametrizes the signal as a Voigtian (con- 
volution of a Breit-Wigner and a Gaussian) 
with mass M w = 0.78256 GeV/c 2 , natural width 
r w = 0.00844 GeV/c 2 and resolution a. The parameter 
s sets the overall strength of the signal. The background 
in each small phase space region was parametrized as a 
fourth order polynomial, 

B(M 3 „, a) = b 4 M^ + b 3 Ml n + b 2 M^ + b x M 3 ^ + b . (4) 

The Q-factor for the event was then calculated as 



Qi 



S(Ml n ,ai 



S(M^,&i) + B(M^,&iy 



(5) 



where M 37T is the event's 3-7T invariant mass and dj are 
the estimators for the parameters obtained from the i th 
event's fit. The signal yield could then be obtained in 
any kinematic bin as 



N 



(0) 



where N is the number of events in the bin. 

The full covariance matrix obtained from each fit was 
used to obtain the uncertainty in Q, <jq. This varied de- 
pending on kinematics; however, it was typically about 
3%. The uncertainty of the extracted yield, in any kine- 
matic bin, was obtained by adding the Q-factor uncer- 
tainties (assuming 100% correlation) to the statistical un- 
certainty of the yield: 



N 

E 



Qi 



' N 

E 



(7) 



Studies were performed using various background 
parametrizations, including polynomials of different or- 
ders, all of which yielded results within the values ob- 
tained for <7q. Therefore, we conclude that no additional 
systematic uncertainty is required. 

Fig. [3] demonstrates the effectiveness of applying this 
procedure in a single center-of-mass energy bin. Fig.[3ja) 
shows the M 37r distribution, integrated over all kinemat- 
ics, and the estimated background using the procedure 
described above. The results are quite plausible; how- 
ever, u) photoproduction provides us with a more strin- 
gent test of this procedure. 
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The distribution of the decay quantity A, which can be 
written in terms of the pion momenta in the u> rest frame 
as 



A = 



\p„+ x p„- 



MAX (\p n + x p„- 1 2 ) ' 



(8) 



must be linear in shape and intersect at A = for u> 



this follows directly from the uj — > n + 7r 7r u 
Fig. ^b) shows the A 



(10) 



events 

amplitude defined in Eq 

distribution, integrated over all kinematics, for events in 
the same bin shown in Fig.[3ja) in the region ±25 MeV/ c 2 
around the to peak, along with the extracted signal and 
background distributions. The signal distribution is well 
described by the function aX. The small discrepancy near 
A = is the result of detector resolution. 

The method we have employed has cleanly separated 
signal from background in the quantity A, even though 
the known linear behavior of the signal was not enforced 
in the fits. In fact, this method has effectively separated 
signal from background in all distributions, successfully 
preserving all kinematic correlations. 

A detailed study of the systematic biases of the back- 
ground subtraction technique was carried out as part of 
this analysis. Not only was the function that was used 
to parametrize the background varied, but the number 
of nearest neighbor events was varied over a wide range 
and several different metrics were used to determine the 
nearest neighbors events. The observed physical mea- 
surements were found to be completely insensitive to 
changes in these parameters over any reasonable set of 
values. Because of this, we associate no additional sys- 
tematic error with these choices. A detailed description 
of this study is contained in Ref. [55] . 

For energy bins near threshold and for "edge" regions 
(i.e. forward- and backward- most angles) in some en- 
ergy bins, the lack of events on both sides of the peak 
leaves the fits unconstrained. In these regions, the en- 
ergy dependence of the Q-factors obtained in the closest 
energy bins for which fitting could be used were projected 
down to the regions in question in order to obtain the Q- 
factors. Fig. [4] shows the results of this procedure in the 
W = 1.735 GeV bin. By studying the A distributions 
in these bins, the systematic uncertainty associated with 
the projected Q-factors in the edge regions is estimated 
to be 5%. In the first two energy bins above threshold, 
the uncertainties are estimated to be 15% and 10% for 
the W = 1.725 GeV and 1.735 GeV bins, respectively. 




M 37I (GeV/c 2 ) 



FIG. 4: The 37r invariant mass distribution in the 
W = 1.735 GeV bin, integrated over all kinematics, for all 
events (unshaded) and for events weighted by the background 
factors, 1 — Q (shaded). 



events were then processed using the same reconstruction 
software as the data. In order to account for the event 
trigger used in this experiment (see Sec. [TTJ) , a study was 
performed to obtain the probability of a track satisfying 
the sector-based coincidences required by the trigger as a 
function of kinematics and struck detector elements. The 
average effect of this correction in our analysis, which re- 
quires three detected particles, is about 5%-6%. 

An additional momentum smearing algorithm was ap- 
plied in order to better match the resolution of the Monte 
Carlo to that of the data. Its effects were studied us- 
ing four-constraint kinematic fits performed on simulated 
"/p — ► ptt + it~ events. After applying the momentum 
smearing algorithm, the same covariance matrix used for 
the data also produced flat confidence level distributions 
in all kinematic regions for the Monte Carlo data as well. 
The simulated to events were then processed with the 
same analysis software as the data, including the one- 
constraint kinematic fits. At this stage, all detector and 
software efficiencies were accounted for. 



V. ACCEPTANCE 

The efficiency of the detector was modeled using the 
standard CLAS GEANT-based simulation package and 
the Monte Carlo technique. A total of 200 million events 
was generated pseudo-randomly, sampled from a phase 
space distribution. Each particle was propagated from 
the event vertex through the CLAS resulting in a simu- 
lated set of detector signals for each track. The simulated 



In order to evaluate the CLAS acceptance for the 
"fp — ► pcu reaction, all kinematic correlations between the 
final state particles must be accurately reproduced by 
the simulated data. Typically, this is done by using a 
physics model when generating Monte Carlo events. Due 
to the lack of any pre-existing precise polarization mea- 
surements in the kinematic regions that contain most of 
our data, this was not an option. Instead, we chose to ex- 
pand the scattering amplitude, M, in a very large basis 
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of s-channel waves as follows: 



M 



mi ,m~ ,m 



~ ^ ^ ^ ^ ^m i7 m-y ,mf tt ) J (^) 

j=ip=± 

where d? denotes a vector of 108 fit parameters, x de- 
notes the complete set of kinematic variables describing 
the reaction, rrii , m 7 , m f , m u are the spin projections on 
the incident photon direction in the center-of-mass frame, 
and A are the s-channel partial wave amplitudes. 



The 



"7r amplitude, which is included in the 



^4's above, can be written in terms of the isovectors and 
the 4-momenta of the pions, /„• and p v respectively, as 
well as the w 4-momentum {q) and polarization (e) as 



Xe^appZ+pZ-P^iq,™,^), (10) 

which is fully symmetric under interchange of the three 
pions. For this reaction, where all final states contain 
oj — > 7T + 7r _ 7r , the isovector triple product simply con- 
tributes a factor to the global phase of all amplitudes. In 
the ui rest frame, Eq. (10) simplifies to 



(Ptt+ x P-k-) ■ e(m u ). 



(11) 



The remaining s-channel structure of the amplitudes A, 
as well as the details concerning the fit parameters, is 
described in [21] . 

Unbinned maximum likelihood fits were performed in 
each W bin in order to obtain the estimators a for the 
parameters a in Eq. ([9]). The results of these fits were 
used to obtain a weight, Ii, for each Monte Carlo event 
according to 



E 



(12) 
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FIG. 5: (Color Online) Example fit result in the 
W = 2.005 GeV bin for data (black squares), phase space ac- 
cepted Monte Carlo events (blue dashed line) and phase space 
accepted Monte Carlo events weighted by the fits discussed in 
Sec. [V] (red solid line). The weighted Monte Carlo provides 
an excellent description of the data. 



where N acc (N gen ) is the number of accepted (generated) 
Monte Carlo events in the bin and the I's are the event 
weights discussed above. An accurate physics genera- 
tor would use the factors of / during the event gener- 
ation stage, rather than weighting the accepted events. 
The resulting acceptance calculation would be the same, 
modulo statistical fluctuations. 

The statistical uncertainties in the acceptance calcula- 
tions due to the finite number of Monte Carlo events gen- 
erated in each kinematic bin were obtained from the stan- 
dard binomial distribution. The systematic uncertainty 
in the acceptance calculation is discussed in Sec. |VII| 



VI. NORMALIZATION 



where Xi is the complete set of kinematic variables of the 
i th event. The weighted accepted Monte Carlo events 
fully reproduce the real data in all distributions, includ- 
ing all correlations (see Figs. [5] and |6|. We note here that 
the results of these fits are not interpreted as physics, 
i.e. they are not considered evidence of resonance con- 
tributions to the uj photoproduction reaction. They are 
simply used in order to provide a complete description of 
the data. 

The acceptance in any kinematic bin is then obtained 

as 

N acc 

E ii 

acc ( f ) = Jf » (I 3 ) 

ff I j 



The measured rate of electrons detected by the tag- 
ger was used to compute the number of photons incident 
on the target by sampling tagger hits not in coincidence 
with the CLAS. These rates were integrated over the live- 
time of the experiment in order to obtain the total pho- 
ton flux associated with each tagger element. Losses of 
photons in the beam line due to effects such as beam 
collimation were determined during dedicated runs us- 
ing a total-absorption counter placed downstream of the 
CLAS |32| . 

The standard electronics hardware that has tradition- 
ally been used to determine the experimental live-time 
malfunctioned during the "glla" data taking period. A 
downstream device used to measure the beam current 
during electron runs |19j was used instead. The relatively 
low count rate in this device during photon running led 
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FIG. 6: (Color Online) <f) versus cos# in the Adair frame (see Sec. VIII B) in cos^",^ bins: Example fit result in the 
W — 2.005 GeV bin for data (top row), phase space accepted Monte Carlo events (middle row) and phase space accepted 
Monte Carlo events weighted by the fits discussed in Sec. [V] (bottom row). Panel (a) shows backwards angles, while panel 
(b) shows forward angles. Note that the weighted Monte Carlo provides an accurate description of the data (x 2 /degrees-of- 
freedom ~ 0.6). 
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to increased uncertainty in the live-time measurement. 
The stability of normalized u> yields for runs with differ- 
ent beam currents was used to estimate this uncertainty 
to be about 3%. 

As was stated in Sec. [TTJ only 40 of the 61 timing el- 
ements of the photon tagger were included in the event 
trigger. Events associated with hits in the "untriggered" 
counters, 41-61, were only recorded if a random hit in 
counters 1-40 occurred during the trigger time window. 
The electron rates used to measure the photon flux, dis- 
cussed above, were used to calculate the probability of 
such an occurrence, Ptrig — 0.467. The measured flux 



for tagger counters 41-61 was scaled down by Pi 



trig 



to 



account for the event trigger. 

Defective electronics in one of the tagger channels led 
to inaccurate flux measurements in the energy bins at 
W = 2.735 and 2.745 GeV. The flux in the energy bin 
aXW= 1.955 GeV was also deemed unreliable due to its 
inclusion of events associated with both triggered and un- 
triggered tagger counters. Differential cross sections are 
not reported in these three energy bins; however, spin 
density matrix elements, which do not require normal- 
ization information, are reported. 



VII. SYSTEMATIC UNCERTAINTIES 

The uj photoproduction cross section, for the case with 
an unpolarized beam and an unpolarized target, must be 
isotropic in the azimuthal angle. Thus, the acceptance- 
corrected u) yields, each obtained in an individual CLAS 
sector, must be consistent with each other. By exam- 
ining the consistency of these yields, we estimated the 
relative uncertainty in the acceptance correction to be be- 
tween 4%-6%, depending on center-of-mass energy. This 
is added in quadrature with uncertainties due to parti- 
cle identification (1.3%) and confidence level (3%) cuts 
to obtain an overall estimated acceptance uncertainty of 
5%-7%. 

It is common practice in photoproduction experiments 
to check the quality of the normalization calculation by 
computing the single pion cross section and comparing it 
to the world's data; however, the two-track trigger used 
in this experiment does not permit such a calculation. 
In order to check our normalization, cross sections were 
also computed for several other reactions from the "glla" 
data set and compared to previously published CLAS 
data. The run-to-run consistency of the normalized u> 
yield was also examined. Based on these studies, we es- 
timate the normalization uncertainty to be 7.3%. When 
combined with contributions from photon transmission 
efficiency (0.5%) and live-time (3%), the total estimated 
normalization uncertainty is 7.9%. 

The acceptance and normalization uncertainties dis- 
cussed above were then combined with contributions 
from target density and length (0.2%), along with 
branching fraction (0.7%) to obtain a total uncertainty, 
excluding contributions from signal-background separa- 



tion that are calculated "point-to-point" , of about 9%- 
11%. In the first two energy bins above threshold, the 
additional uncertainties in the signal-background separa- 
tion method (see Sec. IV I increase this number to 13%- 
17%. 



VIII. RESULTS 



Differential cross sections 



Differential cross sections, da/d cos 0£ m _, were com- 
puted in 109 10-MeV wide bins in W. Each energy bin 
was divided into 20 bins in cos0£ m of width 0.1, al- 
though results could not be extracted in every bin due 
to limitations in the detector acceptance. In total, 1960 
cross-section points are reported here. The centroid of 
each bin is reported as the mean of the range of the 
bin with nonzero acceptance. The results are shown in 
Figs. [7j [H [9] an d [10] The error bars contain the un 



certainties of the yield extraction, discussed in Sec. IV 



along with statistical uncertainties from the Monte Carlo 
acceptance calculations. The overall systematic uncer- 
tainty, discussed in Sec. |VIIt is estimated to be between 
9%-ll%, depending on center-of-mass energy. 

In the "transverse direction" , which we can loosely de- 
fine as | cos0^ m | < 0.8, there are several prominent fea- 
tures present in the data. Near threshold, the transverse 
cross section is mostly flat. Around W ~ 1.9 GeV it 
begins to develop a humped shape and by W ~ 2.1 GeV 
the cross section has two dips. In a concurrent article, 
we present partial wave analysis results obtained from 
this data which attribute these features to various baryon 
resonance contributions [16 . For now, we simply aim 
to draw attention to some of the prominent structures 
present in our measurements. 

A very prominent forward peak begins to rise just 
above threshold and continues to be the dominant fea- 
ture of the cross section up through our highest energies. 
This type of behavior typically indicates the presence 
of strong f-channel contributions. Models of u> photo- 
production, e.g. Refs. [5J 03 0, typically associate 
this peak with i-channcl contributions from 7r°, rj and 
Pomeron exchange. A backwards peak begins to emerge 
around W ~ 2.2 GeV, whose prominence increases as the 
energy increases (although it is always at least one order 
of magnitude smaller than the forward peak) . This could 
be indicative of the presence of contributions in the it- 
channel. Many models of this reaction attribute this peak 
to u-channel nucleon exchange [5J [7J |5J [3] ; however, com- 
parisons of the spin density matrix elements predicted 
by these models to the new high-precision measurements 
presented in this article casts doubt on the validity of 
these models (see Sec. VIII C I. 
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FIG. 7: d co ^gu, — (t^b) versus coa0" m .: Differential cross-section results for bins in the energy range 1.72 GeV < W < 2.00 GeV. 
The centroid of each 10-MeV wide bin is printed on the plot. The lack of reported data points in the W = 1.955 GeV bin is 
discussed in Sec. |VI| The error bars, which do not include systematic uncertainties, are discussed in the text. The additional 
near-threshold background separation uncertainties, discussed in Sec. |IV[ are clearly visible in the first four center-of-mass 
energy bins. Note that the vertical scales are linear up to W of 1.945 GeV and logarithmic above that. 
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B. Spin density matrix elements 

The polarization of the u> can be studied by examining 
the distributions of its decay products. Since the u> is a 
spin-1 particle, its spin density matrix has nine complex 
elements; however, parity, hermiticity and normalization 
reduce the number of independent elements (for an un- 
polarized beam) to four real quantities (of which, three 
are measurable). Traditionally, these are chosen to be 
Pq , p°_ 1 and Re(pi Q ). Our results cover a large range 
of energies and angles; thus, we chose the quantization 
axis to be the photon direction in the overall cm. frame, 
known as the Adair frame [M] . 

The spin density matrix elements can be written in 
terms of the production amplitudes A (i.e. the scattering 
amplitudes M. introduced in Sec. |V| without the lo decay 
piece), as 



P°MM' 



1 E 

m 7 ,m, .rrtf 



A 



mi .m 



.M 



A 



.M' i 



(14) 



where the M, M' refer to the spin projection of the to (on 
the photon direction in the cm. frame) and 



N : 



E El-* 



,M 



(15) 



A I 



is a normalization factor. Using the production ampli- 
tudes obtained from the event-based fits described in 
Sec. [Vj the spin density matrix elements were projected 
out of the partial wave expansion at 2015 (W, cos#" m ) 
points. These data points correspond to the centroids 
of the bins for which cross-section results are reported, 
along with additional points in the W — 1.955 GeV, 
2.735 GeV and 2.745 GeV center-of-mass energy bins for 
which cross-sections results are not reported due to nor- 



malization issues (see Sec. VI I. 

Figs. 11 12 13 and 14 show the p a MM i results extracted 
using the partial wave expansion technique. The error 
bars are purely statistical. The spin density matrix el- 
ements do not rely on normalization information; thus, 
only the acceptance can contribute to the systematic un- 
certainty. Possible effects due to systematic problems 
in the acceptance calculation were examined by analyz- 
ing decay distributions distorted by our estimated accep- 
tance uncertainty. Based on this study, we estimate the 
systematic uncertainties in our results to be as follows: 



cr 00 = 0.0175 
<j x _i = 0.0125 
(710 = 0.01. 



(16a) 
(16b) 
(16c) 



Over most of our kinematics, these results represent 
the first high-precision measurements of p° MM i for lo pho- 
toproduction. Near threshold and at forward angles, the 
cross section develops a strong forward peak, which is in- 
dicative of ^-channel contributions. In this same region, 
the diagonal p^ element decreases sharply as the energy 



increases, or equivalently, as the forward peak increases 
in significance. This is typical of exchange of a spin-0 par- 
ticle in the i-channel where the lo is forced to carry the 
spin of the photon at forward angles. This new precise 
polarization information should help determine the rela- 
tive strengths of the scalar and pseudoscalar exchanges 
(see Sec |yTnC] ). 

At higher energies, starting near W ~ 2.1 GeV, a dip 
in Pq appears at cos6>" m ~ 0.4, which continues to in- 
crease in prominence until about W ~ 2.5 GeV. Above 
this energy, its significance slowly decreases; however, it 
is still present at our highest energies. This dip is located 
near where the forward peak (typically associated with t- 
channcl contributions) has decreased in significance such 
that it is approximately the same size as the cross sec- 
tion in the region < cos0" TO < 0.4. Thus, it is pos- 
sible that this dip results from interference between the 
t-channel and larger-angle production mechanisms. In 
the kinematic regions where the cross section possesses 
the humps and dips discussed in Sec. |VIIIA"| there are a 
number of interesting features found in the spin density 
matrix elements as well. The partial wave analysis we 
performed on this data found that these features are well 
described by baryon resonance contributions |16j . 



C. Interpretation of the data 

In the low-energy regime, these new measurements 
have been used to carry out a mass-independent par- 
tial wave analysis of the reaction jp — > top. The results 
of this analysis, which are presented in a concurrent arti- 
cle |16j and are not discussed in detail here, show clear ev- 
idence of s-channel resonance contributions. This PWA, 
the results of which are different from previous analy- 
ses |HJ Unj El El H2 H3, was the first to benefit from 
the strong additional constraints provided by the high- 
precision polarization results obtained from these data. 

The high-energy measurements have been compared 
to two existing models for lo photoproduction. The 
first is the model of Oh et al |5] which includes pseu- 
doscalar meson (ir° and 77) and Pomeron exchange in the 
t-channel, along with nucleon exchange in both the s- 
and u-channel. It also includes s-channel contributions, 
which are necessary to describe the data in the central 
region of the angular range. The second model is that 
of Laget H| which includes t- and u-channel contri- 
butions similar to that of Ref. [8], but also allows for 
a contribution from two-gluon exchange. In this latter 
model, the two-gluon term is required to describe the <f> 
photoproduction data. 



Fig. 15 shows comparisons of these models to our data 
at W = 2.8 GeV. Both models do a reasonable job 
of reproducing the cross-section measurements. The t- 
channel terms drive the very forward-angle data where 
the agreement is very good. At backwards angles, where 
the it-channel terms dominate, the agreement is not as 
good as it is at forward angles. In the central region, 
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both models agree with the overall shape of the cross 
section; however, the finer structure in the data is not 
reproduced. 

Neither model is able to reproduce the new high- 
precision spin density matrix element measurements pre- 
sented in this article. While some regions are reason- 
ably well described by one model or the other, neither 
gives anything close to good overall agreement. Perhaps 
the most striking discrepancy is that at forward angles, 
where the cross sections are described very well by both 
models, neither provides an excellent description of the 
spin density matrix elements. The high-precision mea- 
surements presented in this article clearly provide new 
stringent constraints on both the nature of the produc- 
tion mechanisms in the high-energy regime, as well as on 
the search for missing baryon resonances. 



D. Comparison to previous measurements 

Previous experimental measurements that overlap our 
energy range have been made at CLAS in 2003 [I], 
at SAPHIR in 2003 0, at Daresbury in 1984 [3] and 
1977 0, and at SLAC in 1973 pQ. Below we compare 
our measurements with each of these previous results. 
The cross sections will be examined first, followed by the 
spin density matrix elements. 

Fig. [16] shows a comparison of the cross-section re- 
sults presented in this article with previously published 
results from CLAS gj and Daresbury [1 |3]. The 
previous CLAS results, four energy bins in the range 
2.624 GeV < W < 2.87 GeV, cover virtually the same 
angular range as the current results. The agreement is 
very good for cos#" m > —0.1; however, there is a siz- 
able discrepancy in the backward direction. At the time 
of the earlier CLAS measurement, the u> polarization had 



only been measured in the forward direction (see Fig. 19 1; 
thus, these values of the spin density matrix elements 
were used in the acceptance calculation. Our results show 
that the polarization is quite different at backward and 
forward angles. Near the edges of the CLAS acceptance, 
e.g. in the backward direction, an incorrect description of 
the polarization can lead to large errors in the acceptance 
calculation. This is most likely the cause of the discrep- 
ancy in the cross sections. The Daresbury results, which 
were only published in the very forward and backward 
regions, are in good agreement with our measurements. 

For W < 2.4 GeV, the previous large acceptance re- 
sults come from SAPHIR [5 . Fig. 17 shows a compari- 
son of the SAPHIR cross-section measurements with the 
present CLAS results. The error bars shown for the 
SAPHIR points do not include systematic uncertainties. 
The agreement is fair, but there are some discrepancies. 
The SAPHIR experiment had better angular coverage; 
however, the CLAS results are more precise. In the for- 
ward direction, the agreement is very good at all energies. 
At moderate angles, \cos8% m \ < 0.5, the agreement is 
good at lower energies but the CLAS results tend to be 



higher as the energy increases. In the backward direc- 
tion, where the CLAS has acceptance, the CLAS points 
are almost always higher than the SAPHIR points. 

Previous spin density matrix element measurements 
are much less precise. The only results published for 
W < 2.4 GeV come from SAPHIR 0. Fig. [18] shows 
a comparison of the SAPHIR results, which consist only 
of four energy bins, each with two angular points, and 
the present CLAS results, which include 1371 total data 
points in this energy range. We note here that the 
SAPHIR collaboration published their results in both the 
Gottfried- Jackson and Hclicity frames, with each mea- 
surement constituting an independent fit to their data. 
Both results were rotated into the Adair frame for com- 
parison. Overall, the SAPHIR results are in good agree- 
ment with our measurements. 

At higher energies, previously published results only 
exist at very forward angles. Fig. [19] shows a compari- 
son of our forward high energy results with those from 
Daresbury [5] and SLAC pQ. The agreement is good. 
For W > 2.4 GeV, the results presented in this ar- 



cos ffi 



Pmm' 
< 0.8. 



IX. CONCLUSIONS 

In summary, experimental results for uj photoproduc- 
tion from the proton have been presented in the energy 
regime from threshold up to W = 2.84 GeV. Both differ- 
ential cross section and spin density matrix element mea- 
surements are reported. The cross-section results are the 
most precise to date and provide the largest energy and 
angular coverage. The results are in fair to good agree- 
ment with previous experiments. For W < 2.4 GeV, 
we present 1181 p Q MM , data points; the previous world's 
data consisted of 8 points. At higher energies, we have 
made the first spin density matrix element measurements 
for cos#" TO < 0.8. Our p° MM i measurements are in 
good agreement with the, rather sparse, existing data. 
The 1960 (W,cos9% m ) cross-section points, along with 
the 2015 (W, cos0"m .) s P m density matrix element data 
points can be obtained at Ref. [25 . 

These new data will have a large impact on our cur- 
rent understanding of vector-meson photoproduction, as 
well as provide a crucial data set in the search for miss- 
ing baryon resonances. A mass-independent partial wave 
analysis performed on these data, which is the first such 
analysis to benefit from the strong constraints provided 
by high-precision polarization information, found strong 
evidence for baryon resonance contributions [16:. Fur- 
thermore, none of the current models of high-energy uj 
photoproduction are able to describe the precise spin 
density matrix element measurements presented in this 
article. We look forward to seeing what impact these 
new results will have on future models of vector-meson 
photoproduction . 
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The centroid of each 10-MeV wide bin is printed on the plot. The error bars, which do not include systematic uncertainties, 
are discussed in the text. 
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FIG. 10: 



(fj,b) versus cos Qt, m , : Differential cross-section results for bins in the energy range 2.56 GeV < W < 2.84 GeV. 



The centroid of each 10-MeV wide bin is printed on the plot. The lack of reported data points in the W = 2.735 GeV and 
W = 2.745 GeV bins is discussed in Sec. |VI| The error bars, which do not include systematic uncertainties, are discussed in 
the text. 
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FIG. 11: (Color Online) (P MM i versus cos#" m : Spin density matrix element measurements, in the Adair frame, for bins in 
range 1.72 GeV < W < 2.00 GeV. The black squares are pot), the red circles are and the blue crosses are Re(p%). 
centroid of each 10-MeV wide bin is printed on the plot. The error bars do not include systematic uncertainties. 
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FIG. 12: (Color Online) p° MM i versus cos#" m : Spin density matrix element measurements, in the Adair frame, for bins in 
range 2.00 GeV < W < 2.28 GeV. The black squares are poo> the re d circles are p?_i and the blue crosses are Re(p%). 
centroid of each 10-MeV wide bin is printed on the plot. The error bars do not include systematic uncertainties. 
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FIG. 13: (Color Online) (P MM i versus cos#" m : Spin density matrix element measurements, in the Adair frame, for bins in 
range 2.28 GeV < W < 2.56 GeV. The black squares are pot), the red circles are and the blue crosses are Re(p%). 
centroid of each 10-MeV wide bin is printed on the plot. The error bars do not include systematic uncertainties. 
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FIG. 14: (Color Online) p° MM i versus cos#" m : Spin density matrix element measurements, in the Adair frame, for bins in 
range 2.56 GeV < W < 2.84GeV. The black squares are poo, the red circles are p?_i and the blue crosses are Re(p%). 
centroid of each 10-MeV wide bin is printed on the plot. The error bars do not include systematic uncertainties. 
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FIG. 15: A comparison of the theoretical models of Oh et 
al [H] (red dashed lines) and Laget [HI [7] (blue dot-dashed 
lines) to the W = 2.805 GeV data presented in this article 
(black squares). 
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FIG. 16: (Color Online) dc J%„ (jib) versus cos6>^. m .: 
Comparison of the present CLAS results (blue circles) with 
previously published results from CLAS !J. (black open 
crosses, triangles, squares and circles) and Daresbury 2, 
[3] (red open diamonds and stars) in the energy range 
2.62 GeV < W < 2.96 GeV. The Daresbury (1977) points 
have no error bars; the points were extracted from a portable 
document format (PDF) image. 
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FIG. 17: (Color Online) d cos d ^, — ()J.b) versus cos8^ m : Comparison of the present CLAS results (blue circles) with previously 
published results from SAPHIR 5. (black open squares) . 
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FIG. 18: (Color Online) p\ [M i in the Adair frame versus cos#" m .: Comparison of the present CLAS results (blue circles) 
with previously published SAPHIR [5] results (open red squares and triangles). SAPHIR extracted results independently in 
the Gottfried-Jackson and Helicity frames — both presented here rotated to the Adair frame. 
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FIG. 19: (Color Online) 

Pmm' i n t ne Adair frame versus coa8^, n : Comparison of the present CLAS results (blue circles) with 
previously published Daresbury [3] (open red circles and triangles) and SLAC [I] (open black squares). 



